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Abstract. The possibility of a correlation between the 
radio (cm)- and 7-ray luminosity of variable AGN seen 
by EGRET is investigated. We performed Monte-Carlo 
simulations of typical data sets and applied different cor- 
relation techniques (partial correlation analysis, x^-test 
applied on flux-flux relations) in view of a truncation bias 
caused by sensitivity limits of the surveys. For K-corrected 
flux densities, we find that with the least squares method 
only a linear correlation can be recovered. Partial corre- 
lation analysis on the other side provides a robust tool to 
detect correlations even in flux-limited samples if intrinsic 
scatter does not exceed ~ 40 % of the original 7-ray lumi- 
nosity. The analysis presented in this paper takes into ac- 
count redshift bias and truncation effects simultaneously 
which was never considered in earlier papers. 

Applying this analysis to simultaneously observed 
radio- and 7-ray data, no correlation is found. However, 
an artificial correlation appears when using the mean flux. 
This is probably due to the reduction of the dynamical 
range in the flux-flux relation. Furthermore, we show that 
comparing the emission in both spectral bands at a high 
activity state leads to no convincing correlation. 

In conclusion, we can not confirm a correlation be- 
tween radio and 7-ray luminosities of AGN which is 
claimed in previous works. 

Key words: galaxies: active - quasars: general - gamma 
rays: observations - Methods: statistical 
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1. Introduction 

Many attempts have been made in the past to investigate 
correlations between radio (cm)- and 7-ray luminosities of 
AGN (Sleeker et al. 1993, Padovani et al. 1993, Salamon 
& Sleeker 1994). By applying regression and correlation 
analysis to the 2.7 GHz- and 5 GHz- luminosity and the 
7-ray emission in the EGRET energy band, a (nearly) lin- 
ear correlation was found. The use of luminosities instead 
of fluxes, however, always introduces a redshift bias to 
the data, since luminosities are strongly correlated with 
redshift. Especially in samples which cover a wide range 
of distance, a correlation will appear in a luminosity plot 
even when there is no correlation in the corresponding flux 
densities (Elvis et al. 1978). Feigelson & Berg (1983) show 
that, if there is no intrinsic luminosity- luminosity corre- 
lation, no correlation will appear in the flux-flux relation. 
On the other hand, though the redshift dependence can 
be removed, intrinsic correlations between luminosites Li 
and L2 may be lost in the flux diagrams /1-/2: if Li ^ L^, 
then /i - fi^z^^^"^\ For B 7^ 1 each value wiU be 
multiplied by a 'random' value z'^^^~^^ causing any intrin- 
sic Li — L2 correlation to be smeared out. The apparent 
correlation is maintained only when the underlying rela- 
tionship is linear (Feigelson & Berg 1983). Furthermore, 
any observational uncertainties will hamper the finding 
of correlations using fiux diagrams. It is therefore crucial 
to estimate the influence of the redshift bias on the cor- 
relations between the emission from the two wavebands. 
Non-parametric partial correlation coefficients (Kendall's 
T or the Spearman rank correlation coefficient Rs) have 
been used to deal with this problem (Dondi & Ghisellini 
1995). One limitation of these rank correlation tests is the 
failure to take into account possible changes in rank due 
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to observational uncertainties. 

Observational flux-limits of the samples arc another seri- 
ous bias which influences any correlation analysis as they 
restrict the populated region in the luminosity-lmninosity 
diagram to a narrow band. Therefore, Fcigclson & Berg 
(1983) proposed to include all upper limits to avoid artifi- 
cial correlations and incorrect conclusions (Schmitt 1985), 
and suggest the use of survival analysis. However, if the 
censored data points arc not distributed randomly, but 
localize a particular area, survival analysis may give mis- 
leading results (Isobe 1989). Furthermore, this analysis 
can not account for a bias caused by misidentifled sources 
or by objects which are completely lost due to the low sen- 
sitivity of the instrument. Those truncation effects must 
be seriously considered when EGRET data are used. 
Another problem arises when considering the nature of 
the sources. Blazars are known to be strongly variable in 
the 7-ray as well as in the radio band on time scales of 
days to months (von Montigny et al. 1995). Therefore, si- 
multaneous observations should be the adequate data for 
a correlation analysis. However, due to the lack of such 
data, the mean (Padovani et al. 1993) or the brightest flux 
values (Dondi & Ghisellini 1995) of the sources have been 
used in the past. For highly variable sources this choice 
reduces the dynamical range in the luminosity-luminosity 
plot essentially, and can mimic a correlation. 

In this paper, we evaluate the reliability of the partial 
correlation analysis in luminosity-luminosity plots as well 
as X^'fits in flux-flux diagrams for flux-limited truncated 
samples with strongly variable sources. For this purpose 
we perform Monte-Carlo simulations of typical data sets 
with known degrees of correlation and flux sensitivities, 
which is described in section 2. No censored data are used. 
The simulations provide a useful test of how correlation 
analysis methods deal with such selection effects. 

The application of a database of EGRET-blazars in 
the radio (cm)- and 7-ray regime is presented in the last 
section. 

Throughout this paper the Hubble constant Hq = 75 
km Mpc~^ and the deceleration parameter qo = 0.5 
have been used. 

2. Investigation of correlation techniques 

2.1. Monte- Carlo simulations and correlation analysis 

We consider a sample of objects with a wide range of dis- 
tances observed in two spectral bands. To approximate 
the real case, we assume a sample of flat-spectrum ra- 
dio quasars, with strong emission in the radio and 7-ray 
regime. The range of the radio luminosity Lr is chosen 
so that the simulated data resemble our database, and is 
set at 40 < logL^f erg s~^] < 46. The radio luminosity 
function Pradio is taken from Dunlop & Peacock (1990) 



Pradio ~ [( 



Lc{z) 



Lr 
Lc{z) 



1.961 



(1) 



where Lc(z) is the 'evolving' break luminosity and was 
parametrized as \ogLc{z) = 25.26 -|- I.I82; — 0.28^;^ with 
the redshift z. 

For the case of a i-, — Lr correlation the 7-ray lumi- 
nosity is calculated from the relation 



logL^ = A + B\ogLR + e{(j) 



(2) 



where A and B are constants and free parameters. The 
term e(cr) is a random noise component following a normal 
distribution with dispersion a. 

If the radio and 7-ray luminosities are assumed to be 
uncorrelated, the probability of detecting an object with 
a luminosity L in the 7-ray regime is calculated from the 
luminosity function p-y of EGRET-blazars (Chiang et al. 
1995) 



P7 
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where 71 = 2.9, 72 = 2.6, Lr = iC^erg s~^/i~^ and h = 
Hq/IOO km s^^ Mpc^^ . A pure luminosity evolution is 
incorporated with L~,,{z) ^ (1 + z)^-^ (sec Chiang et al. 
1995). The 7-ray luminosity range is again taken from the 
observation: 44 < logL^[ erg ,s^^] < 50. 

The objects are distributed in three dimensional space 
according to their evolution properties with redshifts z 
lying between 0.001 < z < 2.5. The differential redshift 
distribution dN{z) ~ PradiodV (z) {dV{z) — the comov- 
ing volume element) expected from the radio luminosity 
function is 
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with Slim the sensitivity limit of the instrument and c the 
speed of light. Using the relation 



Lband = 0.43 • IQ^'Sband ' E{1 + z) 

Al + z)-^Jl + z.2 



\a-l 



Ho 



(5) 



with Ho the Hubble constant. Stand meaning either the 
radio flux Sradio or the 7-ray flux S^ , E the mean photon 
energy and L^and the corresponding 'z^L(i/)' luminosity 
(Weedman 1986), we convert the luminosity to the flux. 
The term (1 -|- z)°'~^ is responsible for the K-correction. 
The average energy spectral indices for the respective 
spectral bands {Sband ~ z^~") are aradio = for the radio 
band, and a-y = 1 for the 7-ray regime, where the latter 
also defines the mean photon energy of E =470 MeV. All 
objects which have fluxes lower than a specified limit are 
dropped from the data set. To account for the uneven ex- 
posure coverage of EGRET, we use a sigmoidal probability 
distribution for the detection of an object: 

PiSiim) = 1 - [1 + exp {C{Sum + So))]-' 
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The parameters C and Sq are chosen to be C = 6.1(20.5) 
and 5*0 = — 1.25(— 0.55) for the radio (7-ray) band. This 
means that the crcdibihty limit in the radio band ranges 
from ^ 0.5 ... 2 Jy which can be considered as a result of 
an overlay of several radio catalogs with different cred- 
ibility limits. The parameters corresponding to the 7- 
ray band are estimated from the flux threshold map of 
EGRET made from phase 1 and 2 of the mission. We cre- 
ated simulated data sets for different correlation slopes (B 
= 0.8, 1.0, 1.3 and 1.5), assuming either sensitivity limits 
as described above or perfectly sensitive detectors (i.e. flux 
densities down to 0.1 mJy in the radio band, and 10"^""^ 



photon s 



in the 7-ray regime are observed) . 



in the radio The simulations were done until each sample 
contains N = 25 and 12 objects, which are typical num- 
bers for the simultaneously observed data sets and the 
flaring state samples, respectively. For each parameter set 
the whole procedure was then repeated 100 times. 

We then analysed the data sets using Spearman's rank 
order partial correlation cocfhcient Rg (Macklin 1982). 
The partial correlation describes the relationship between 
two variables when the third variable is held constant. In 
this way the strong redshift dependence of the radio and 7- 
ray luminosity can be eliminated. To date the behaviour 
of Spearman's correlation coefficient in the presence of 
upper limits caused by sensitivity limits of the data sets 
is not known, in contrast to Kendall's r, a further non- 
parametric correlation coefficient. The parametric meth- 
ods are more efficient if one knows the exact distribution 
of the data. However, in astronomy we do not know the 
exact distribution functions and hence the no n- parametric 
methods are usually preferable. There has also been some 
debate whether correlations can be discovered more se- 
cure by comparing fluxes instead of luminosities ( Elvis 
et al. 1978). In order to examine how accurately flux-flux 
diagrams recover quantitative information from the data 
when flux-limits and uncertainties are taken into account, 
we carried out a x^-test of a linear relation y = a + bx 
between the logarithms of the K-corrected flux densities. 
The flux density values were provided with uncertainties 
of 4% and 20% in the radio and 7-ray band, respectively, 
which are typical values for those spectral bands. The re- 
sults are summarized in Table 1. 

2.2. Results and discussion 

The results of the partial correlation analysis (see Table 
1) show that a correlation can be recovered at a signifi- 
cance level of > 99.5% for both complete and flux-limited 
samples provided that the intrinsic luminosity scatter is 
not too large {a < 0.05). Increasing this noise factor up 
to a = 0.5, a correlation can still be detected at an 95 % 
level for a complete data set. The most striking effect of a 
sensitivity limit is a generally lower correlation coefficient 
compared to complete data sets. High redshift objects can 
only be detected at high powers. Hence, sensitivity limits 
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Fig. 1. Typical distribution of a flux-limited sample of ob- 
jects wliicli possesses no correlation between tfie radio and 
7-ray luminosity. The dashed lines follow the sensitivity limits 
of the sample. 
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Fig. 2. Typical distribution of a complete sample of objects 
which possesses no correlation between the radio and 7-ray 
luminosity. 



exclude mainly high-redshift objects with low luminosi- 
ties. This effect is shown in Fig.l and Fig. 2 where sim- 
ulated uncorrelated data are presented for a flux-limited 
and a complete sample of objects. While low-luminosity 

sources are rarely observed in flux- limited data sets, they 
are more common in samples which possess no flux-limits. 
This causes an overestimation of the redshift dependence 
in flux-limited samples. Therefore flux- limits seem to lower 
the degree of correlation. 

Analysing flux-flux diagrams with the x^-tests of a lin- 
ear relation can accoimt for observational uncertainties in 
contrast to a partial correlation analysis. However, only a 
linear intrinsic correlation (L^ ~ Lj-adio) can be recovered 
with this method no matter whether the sample obeys 
a sensitivity limit or not. Except for the case of a lin- 
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Table 1. Results of the correlation analysis of simulated data: partial correlation analysis using Spearman rank order correlation 
coefficient Rs (with the probabilities (*) of erroneously rejecting the null hypothesis (i.e. no correlation), and a least square 
analysis with the minimum x^, the goodness-of-flt probability (**) and the slope b of the fit as indicated above. 100 trials 

wore carried out with each sample containing N = 25 and 12 sources. In case of a correlation, a random noise component with 
a = 0.05 was chosen. The analysis were carried out between the logarithms of the luminosities and K-corrected flux densities, 
respectively. For the x^-Ati the flux density values were provided with uncertainties of 4% and 20% in the radio and 7-ray band, 
respectively. A correlation is said to exist for a chance probability (*) < 5 %, and a significance (**) > 68 % using the partial 
correlation analysis and least square method, respectively. 



ear correlation the derived slopes of the correlations show 
misleading trends and the goodness of the x^-fit drops 
down when sensitivity limits restrict the flux range. The 
low range of apparent brightness, the 'distance random 
noise factor z-^^^~^^', observational uncertainties and any 
intrinsic random noise factor e((T) cause a large scatter 
which prevents a reliable goodness-of-fit probability. In 
addition, we note that any intrinsic correlation ~ 
will be randomized by a factor (1 + 2;)"''-i-S(ar„jiio-i) 
(^ 3 in our case) if the fluxes are not K-corrected. Espe- 
cially for samples which cover a wide range of redshifts 
but a small range of apparent brightness it is necessary 
to use K-corrccted fluxes. Furthermore, when the intrin- 
sic scatter causes a standard deviation from the expected 
luminosity of > 12% {a > 0.05), we even can not find a 
linear correlation. Therefore, an intrinsic correlation could 
easily disappear using this method. In the case of a cor- 
relation with low intrinsic scatter (ct ^ 0.05), the right 
slopes are recovered within the parameter uncertainties. 



3.1. The sample 

Our sample consists of 38 identified extragalactic point 
sources observed by EGRET between April 1991 to 
September 1993 (phase 1-^2). EGRET covers the high en- 
ergy 7-ray range from about 0.03 to 10 GcV with a field of 
view of about one-half steradian. The typical duration of 
an observation is two weeks. Only sources detected < 25° 
off- axis and with a significance > 4(7 at least once during 
the observation run are used in this analysis. The flux val- 
ues used here are the integral flux (> 100 MeV), and are 
already published in Thompson et al. (1995). 

Multifrequency radio observations of FSRQ with the 
100-m-Effelsberg Telescope have been performed in par- 
allel to the CGRO all-sky survey. All flat spectrum, vari- 
able sources stronger than 1 Jy at 6cm wavelength in the 
50" field of view of EGRET have been observed quasi- 
simultaneously at 2.8cm, 6cm and 11cm. Sources, which 
have been detected by EGRET, were subsequently moni- 
tored at shorter time intervals. The observational method 
is described by Reich et al. (1993). 



3. Application to EGRET blazars 
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3.2. Data compilation 

All sources used in this analysis are known to be vari- 
able. For sources of known redshift, we define Lband as 
the 'pL{uy luminosity at the indicated frequency. With 
a flux S and a redshift z, this luminosity L is given by 
equation (5). For five objects with unknown redshift z, we 
set 2: = 1 which is the mean value of our sample. For all 
sources an energy spectral index of aradio = in the radio 
band and a^^ = 1 in the 7-ray regime is assumed which 
are typical values. 

In Fig. 3 the 7-ray flux is plotted against the simul- 
taneously observed 2.7 GHz radio flux. The simultaneous 
observations were done at the Effelsberg telescope. We 
denote a pair of radio and 7-ray observations as 'simul- 
taneous' if their observation dates do not differ by more 
than three weeks. Repeated EGRET observations of an 
AGN are taken as independent sources. The correspond- 
ing luminosity diagram is presented in Fig. 4. 

We also investigated a possible correlation between the 
8 GHz and the 4.8 GHz and the 7-ray luminosities in a 
flaring state. For the compilation of this data set, we have 
taken for each source the highest radio flux found in the 
literature (AUer et al. (1985), Kiihr et al. (1981), White & 
Becker (1982), Wafl & Peacock (1985), Reich et al. (1993), 
Fiedler et al. (1987), Seielstad et al. (1983), Waltman et 
al. 1991, Wright et al. (1990) and the PKSCAT90 and 
the maximal 7-ray flux observed by EGRET. In order to 
ensure that we really get the peak fluxes, we only use 
the historical flux values for > 100 % variable sources 
(= {Smax - Smin)/Smin with Smin the lowest detection 
or upper limit, and Smax the highest detection ever ob- 
served) . 

The average flux values at 4.8 GHz and 8 GHz were 
determined from the historical minimum and maximum 
found in the literature (reference see above) . The mean 7- 
ray flux was calculated similarly from the minimum and 
maximum detection during phase 1 and 2 of the EGRET 
mission. 
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Fig. 3. Simultaneously observed 7-ray flux (> 100 MeV) (25 
detections and 16 upper limits) versus 2.7 GHz-radio flux deu- 
sity. The dashed line follows the sensitivity limits of the sample. 
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Fig. 4. Simultaneously observed 7-ray luminosity versus 2.7 
GHz-radio luminosity. The dashed line follows the sensitivity 
limits of the sample. 



3.3. Results and discussion 

The results of the correlation analysis between different 
kind of samples arc presented in Table 2. We considered 
simultaneous detections (including upper limits UL) at 2.7 
GHz and 10 GHz as well as data sets containing the re- 
spective highest flux of each source in the radio and 7-ray 
regime. To allow comparison to previous investigations, 
the average flux value of each source in the two spectral 
bands is also correlated. Analysing K-corrected flux-flux 
diagrams a positive signal could not be found for any of 
the considered samples. We therefore conclude that either 
an intrinsic scatter has smeared out a possible correlation 
or at least a linear correlation between the radio (cm)- 
and 7-ray emission does not exist. 



In the last chapter we showed that the Spearman's 

partial correlation coefficient can recover intrinsic corre- 
lations more secure. Nevertheless, even this method gives 
no significant indication for a correlation between simul- 
taneous observations. This is also true when upper limits 
are taken into account and Kendall's partial correlation 
coefficient is used. For this method it is possible to state 
the significance level even for censored data (Akritas & 
Siebert 1996), in contrast to Spearman's Rg. 

The marginal correlation fomid between the flaring 4.8 
GHz- and 7-ray observations looses signiflcance when con- 
sidering the 8 GHz - 7-ray high-activity state relation 
which possess a high chance probability of > 30%. The 
flaring state observations show only at one radio wave- 
length a marginal positive result, which we therefore con- 
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data 


N 


corr.coeff. 


prob.* 


corr. ? 




signif.** 


slope b 


corr. ? 


10 GHz - 7-ray obs.(simult.det.) 


25 


0.158 


0.465 


NO 


352.9 


0. 


-0.109 ±0.029 


NO 


10 GHz - 7-ray obs.(simult.det.-l-UL) 


42 


0.079 


0.288 


NO 










2.7 GHz - 7-ray obs.(siimilt.dct.) 


22 


0.046 


0.751 


NO 


298.5 


0. 


-0.031 ± 0.026 


NO 


2.7 GHz - "i-ray obs.(siiiiult.<l('t. + l'L) 


11 


().()().■', 


().12() 


NO 










max. 4.8 GHz - 7-ray obs. 


12 


0.594 


0.053 


marg. 


225.9 


6 • lO^*-' 


0.450 ± 0.029 


NO 


max. 8 GHz - 7-ray obs. 


11 


0.363 


0.314 


NO 


238.0 


0. 


0.477 ± 0.036 


NO 


4.8 GHz - 7-ray obs. (mean values) 


38 


0.347 


0.035 


YES 


159.3 


2 • IQ-^' 


-0.225 ± 0.041 


NO 


8 GHz - 7-ray obs. (mean values) 


28 


0.405 


0.035 


YES 


144.4 


2-10-^^ 


0.316 ± 0.073 


NO 



Table 2. Results of the correlation analysis of N observed data points: partial correlation analysis using the Spearman rank 
order coefficient Rs or Kendall's t for the case of data sets where upper limits are included (with the probabilities (*) of erro- 
neously rejecting the null hypothesis (i.e. no correlation)), and a least square analysis with the minimum x^, the goodness-of-fit 
probability (**) and the slope b of the fitted line. The analysis were carried out between the logarithms of the luminosities and 
K-corrected flux densities, respectively. 



sider as accidental. Hence, a convincing correlation be- 
tween flaring detections does not exist. 

Relating the mean 4.8 GHz- as well as the 8 GHz-radio 
and 7-ray luminosities gives a 2 cr result for a positive cor- 
relation (see Table 2). However, using averaged flux values 
in highly variable sources induce a bias which could mimic 
a correlation. This can be explained as follows: Suppose, a 
number of variable objects with totally uncorrelated fltix 
values are observed. The flux at the high activity states 
of the sources will be found at the upper right side of the 
correlation plot while the low activity states are crowded 
at the lower left side. After averaging the low and high 
activity data the resulting flux values occupy the narrow 
region in the middle of the correlation plot, and will ar- 
range themselves along a straight line in the luminosity- 
luminosity-plot. This causes a bias which is not accounted 
for by the correlation methods. Hence, reducing the dy- 
namical range in the flux-flux relation by the averaging 
procedure may cause an artificial correlation. 

Our results using Monte-Carlo simulations can not 
confirm earlier findings concerning correlations between 
the radio (cm)- and 7-ray emission of blazars (Stecker 
et al. 1993, Padovani et al. 1993, Dondi & Ghisellini 
1995, Salamon & Stecker 1994). Note, that the correlated 
quantities they used were not observed simultaneously al- 
though the sources are known to be highly variable. In- 
stead, either the brightest (Dondi & Ghisellini 1995) or the 
mean fluxes (Stecker et al. 1993, Padovani et al. 1993) were 
taken. However, the use of average fluxes reduces the dy- 
namical range in the flux diagrams and hence may mimic 
a correlation (see also Table 2). We conclude that our 
negative finding of a correlation between simultaneously 
observed sources is likely to be caused by the large scatter 
of the data points. Note also, that truncation effects as 
well as the influence of the redshift bias were not consid- 
ered simultaneously in earlier works. Instead, mostly sur- 
vival analysis was applied to eliminate biases caused by 
the sensitivity limits of the surveys. This method, how- 



ever, was at then not able to also deal with a bias caused 
by the strong redshift-dependence of the correlated quan- 
tities. On the other hand, Dondi & Ghisellini (1995) have 
taken this effect into account by using partial correlation 
analysis, but neglected a bias caused by flux-limits. We 
therefore conclude that both, truncation effects as well as 
redshift bias, seriously influence the previously performed 
correlation analysis of radio- and 7-ray-loud blazars. 



4. Conclusions 

In previous work a correlation between the radio (cm)- and 

7-ray luminosity of blazars seen by EGRET was claimed. 
For the correlation analysis biases caused by the limited 
sensitivity of the instruments and the strong redshift de- 
pendence of the luminosities must be considered. While 
the effects of flux-limits can be taken into account by 
including upper limits, partial correlation analysis deals 
with the redshift bias. In this paper we considered both 
effects simultaneously in contrast to earlier works. For 
this purpose we performed Monte-Carlo simulations of a 
typical sample of variable blazars with known degree of 
correlation and completeness. We then analyzed the data 
sets by partial correlation analysis using the Spearman 
rank order correlation coefficient and a least square fit in 
flux-fiux relations. The simulations reveal how correlation 
analysis deals with truncation biases caused by a limited 
sensitivity of the instrument. While in flux-flux diagrams 
only linear correlations can be recovered (provided that K- 
correctcd fluxes are used) with the least square method, 
the partial correlation analysis gives correct results even 
for flux-limited sample. This is demonstrated for expo- 
nents B = 0.8 ... 1.5 (L-y ~ L^) with random noise fac- 
tors up to an order of magnitude. Furthermore, in data 
sets with flux-limits, the partial correlation analysis gives 
a stronger redshift dependence of the luminosities com- 
pared to complete samples. 
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Applying the x^"™^6thod on flux-flux diagrams we 
found no evidence for a correlation between the simultane- 
ously observed (2.7 GHz and 10 GHz) or flaring radio flux 
and luminosity (4.8 GHz and 8 GHz) and 7-ray luminosi- 
ties above 100 MciV. The same resiilt is obtained rising 
the partial correlation analysis. The apparently positive 
signal between the average radio- and 7-ray luminosities 
is caused by the restriction of the dynamical range in the 
correlation diagrams. 

We conclude that a correlation between the radio (cm)- 
and 7-ray luminosities of AGN can not be claimed. 

Note however, that there is evidence that both radio- 
and 7-ray emissions from blazars are strongly beamed and 
not isotropic (von Montigny et al. 1995). The apparent 
source luminosity is related to the intrinsic luminosity by a 
factor which depends on the bulk Lorenz factor of the out- 
moving radiation source region within the jet and the an- 
gle between the line of sight and the jet axis. Those quan- 
tities are thought to be different for each object and each 
spectral band. Therefore, the beaming effect may smear 
out a possible correlation between the intrinsic luminosi- 
ties, or worse, could cause a misleading correlation slope. 
Correlating luminosities which are corrected for beaming 
is unfortunately to date not possible since the Doppler fac- 
tors in the 7-ray band are not known from observations. 

The fact that we do not find any radio-7-ray corre- 
lation for both simultaneously observed luminosities and 
high activity state observations has several important im- 
plications. 

One may ask whether the lack of a correlation would 
rule out any of the proposed 7-ray emission processes in 

blazars. Especially the (inhomogeneous) synchrotron-self- 
Compton (SSC) models (Jones et al. 1974, Marscher 1980, 
Konigl 1981, Marscher & Gear 1985, Ghisellini & Maraschi 
1989, Marscher & Bloom 1992, Maraschi et al. 1992) seem 
to predict a correlation between the synchrotron emission 
and the inverse Compton component, as the synchrotron 
photons are upscattered to 7-ray energies by the same 
beamed relativistic electrons. However, there are several 
reasons why we do not expect to find any relation between 
the luminosities in both spectral bands, even if we con- 
sider the SSC mechanism as the main 7-ray emission pro- 
cess. First, rapid variability, especially in the 7-ray regime, 
causes a large scatter in the luminosity and flux diagrams, 
and prevents recovering the true correlation slope even 
if a correlation would exist. Second, different relativistic 
Doppler factors D for the different objects induce an addi- 
tional scatter since the observer considers different parts of 
the spectrum depending on the beaming factor D. Third, 
different 7-ray emission components may contribute differ- 
ently for each individual object contaminating the relation 
with the radio band (i.e. the lack of strong emission lines in 
BL Lac objects indicates that the 7-ray emission induced 
by inverse Compton scattering of photons produced in the 
broad line region is small (Sikora et al. 1994)). Note also, 
that the radio emission may be influenced by the con- 



tribution of a quiescent component (Valtaoja et al. 1988). 
Fourth, the different variability behaviour in the radio and 
7-ray band points to a different spatial origin of the two 
luminosities which excludes a direct link. Fifth, there is 
evidence for a time lag between the radio and 7-ray flare 
(Miicke et al. 1996, Valtaoja & Terasranta 1995, Reich et 
al. 1993). Simultaneously observed luminosities arc there- 
fore not expected to obey a direct correlation. The same 
is true for the historical maximum relations between the 
radio and 7-ray regime since those flaring states must not 
be necessarily physically connected. Instead, the fluence 
of the two corresponding radio and 7-ray flares may be 
related. However, a single radio outburst can often not 
be isolated since radio flares often overlap especially at 
low cm-frequencies. Unfortunately, the lack of long-term 
7-ray light curves of blazars does not enable a meaningful 
quantitative study at present. Last, the negative finding 
of a luminosity correlation could still be compatible with 
a similar energy content of the radio and 7-ray producing 
particles. Thus, we do not consider the non-correlation 
between simultaneously observed luminosities as strong 
evidence against the SSC model. 

However, the non-correlation yields other astrophysi- 
cal consequences. For both, the local luminosity density of 
blazars and their evolution properties, there is no hint for 
a direct link between radio (cm)- and 7-ray energies. So, 
the 7-ray luminosity function can not be predicted from a 
known radio luminosity function and its evolution as has 
been used by several authors (Stecker et al. 1993, Salamon 
& Stecker 1994, Padovani et al. 1993, Stecker & Salamon 
1996) in order to estimate the contribution of radio-loud 
AGN to the extragalactic diffuse 7-ray background. In- 
stead, we suggest to model the 7-ray background induced 
by unresolved blazars by relating the energy content of 
particles to the 7-ray emission in AGN and take into ac- 
count cosmological distances and the distribution of beam- 
ing angles, duty factors, etc.. This will be subject of future 
work. 
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